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Abstract: We show by numerical modeling that saturation of the 

population inversion reduces the stimulated thermal Rayleigh gain relative 

f-^ • to the laser gain in large mode area fiber amplifiers. We show how to 

.^ \ exploit this effect to raise mode instability thresholds by a substantial 

factor. We also demonstrate that when suppression of stimulated Brillouin 
scattering and the population saturation effect are both taken into account, 
counter-pumped amplifiers have higher mode instability thresholds than 

Y^ \ co-pumped amplifiers for fully Yb^+ doped cores, and confined doping can 

further raise the thresholds. 
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1. Introduction 

In earlier papers we described a stimulated thermal Rayleigh scattering process (STRS) that can 
account for observed modal instability in large mode area fiber amplifiers, and we described in 
detail our numerical model of this process ([l}j5]. Our model computes the mode coupling gain 
and the laser gain for fiber amplifiers with practically achievable step index profiles operating 
at realistic pump and signal powers. It was shown that the quantum defect heating associated 
with laser amplification, in conjunction with a frequency offset between the strong fundamental 
mode and the weak parasitic higher order mode, leads to a high exponential gain for the parasitic 
mode. Consequently, mode instability has a sharp power threshold, above which the output 
beam quality is severely degraded. 

The essence of the STRS process responsible for mode instability is that laser gain necessar- 
ily deposits quantum defect heat in the core of the amplifier fiber, and the asymmetric heating 
produced by the asymmetric signal irradiance profile due to interference between modes LPqi 
and LPii leads to an asymmetric thermal lens that couples light between those two modes. 
An additional requirement is that there be a phase shift between the temperature grating and 
the signal irradiance grating which, in our STRS model, is created by the time lag between an 
irradiance grating traveling along the fiber and the temperature grating that it creates. 

Our model imposes a steady-periodic condition on the temperature grating because we as- 
sume the frequency offset between modes has a narrow linewidth. Alternative STRS models 
have been presented by Hansen et al. ||6l|2l and by Dong IH. They used similar approxima- 
tions, including the steady-periodic assumption, and they predict instability thresholds similar 
to ours. Another model by Ward et al. Q is based on STRS, but without the steady-periodic 
assumption. It also predicts thresholds similar to ours. 

The primary difference in the physics of the various models is that the models of Hansen et 
al. and Dong assume the profile of the quantum defect heating matches the profile of the signal 
light, while our model and that of Ward et al. compute the heat profile based on the local change 
in either the pump or signal irradiance calculated using the local upper state population. In the 
model of Ward et al. the increase of signal power in a mode due to laser gain is found from 
the overlap of the local gsm g{x^y^z) with the field of that mode. In ourBPM model laser gain 
increases the total signal field locally and is then apportioned among the modes automatically 
by diffraction in the presence of the core index step. Either the local signal field growth or the 
local pump irradiance loss is used to compute the quantum defect heating. 

Hansen et al. and Dong showed that using their assumed heat profile the mode coupling 
gain is related in a simple way to the laser gain. This would imply that the mode instability 
threshold is largely determined by the net laser gain, and can be adjusted only by changing the 
modal profiles and their overlap with the Yb^+ doping profile. In fact, they show that by using 
a low value of the V parameter or by confining the Yb^+ doping to the central portion of the 
core, the threshold is raised. However, we show in this report that the existence of depletion 
of the upper state Yb^+ population breaks the simple connection between laser gain and mode 



coupling gain, making it possible to design fibers with substantially higher mode coupling 
thresholds than predicted by Hansen et al. and by Dong. An amplifier that is designed to reach 
a specific level of laser amplification can be designed to avoid the mode instability. Specifically, 
by designing the amplifier so it has a higher ratio of pump cladding diameter to core diameter, 
the STRS threshold can be raised. Of course this requires a longer fiber and a smaller core, a 
combination that is problematic for SBS suppression. In practice the amplifier must achieve a 
suitable balance between STRS and SBS suppression. We will also discuss how to achieve this 
balance. 

2. Transverse hole burning 

In computing the upper state population fraction «„ we use the steady state expression 



riui^.y) 



Ipa"/hvp + h{x,y)a'^/hvs 



Ip{a" + ai)/hv„ + h{x,y){af + (7/)//iV, + 1/t ' 



(1) 



Here the cj's are the absorption and emission cross sections for the pump and the signal, the 
v's are the optical frequencies, and Ip and Is are the pump and signal irradiances. The amplifier 
parameters used throughout this report are listed in Table 1 . We model a fiber with a 50 /im 
diameter step index core with a typical numerical aperture of 0.054. The Yb^+ doping density 
is also typical of high power amplifiers. 





Table 1 . Amplifier 


parameters 


"core 


50 /im 


^dope 


30-50 jum 


dclad 


100-500 Aim 


Nvh 


S.OxlO^^m-^ 


K 


976 nm 


A. 


1032 nm 


^P 


2.47 X 10-24 m2 


^P 


2.44X 10-24 m2 


q" 


5.80X 10-27 m2 


< 


5.0x 10-25 m2 


Pp 


varies 


Ps 


low 


dn/dT 


1.2x10-5 


L 


varies 


P 


2201 kg/m^ 


C 


702 J/kg-K 


^core 


1.451 


ndad 


1.45 


T 


901 jUs 


K 


1.38W/m-K 


NA 


0.054 


V 


8.2 


Aeff(LPoi) 


1175^m2 







Figures[T]and|2]show upper state population profiles, nu{x,y) computed at three points along 
co-pumped and counter-pumped fiber amplifiers with a cladding diameter of 400 /im. The 
signal light is all in LPqi here. Near the input end of the co-pumped fiber the pump is strong and 
the signal is weak so the upper state population is weakly depleted. The shape of the population 
depletion closely matches the signal irradiance profile here. The undepleted population in the 
region with I^ = is «„ = a"/ [a" + Op) « 0.5. Further along this fiber, near the crossing of the 
signal and pump powers, the saturation is strong at the mode center but moderate near the core 
boundary, while near the output end saturation is strong across the entire core. If the same fiber 
is counter-pumped the signal and pump powers are approximately equal along the full length of 
the fiber, so the saturation resembles that of the middle diagram of the co-pumped fiber where 
the signal and pump powers are equal. The degree of saturation would be reduced in a fiber 
with a smaller cladding, and increased in one with a larger cladding. 

The quantum defect heating is proportional to the pump absorption, which is determined 
from the value of nu{x,y) and the pump irradiance. The latter is assumed uniform across the 
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Fig. 1. Signal and pump powers versus z for co-pumped fiber with dcoie = (^dope = 50 /xm, 
"^clad = 400 lira, operating at the mode instability threshold. The inset figures show the 
upper state fraction n„ profiles at z = 0.1, 1.4, and 4.0 m. The range of the «„ axis is 0-0.5. 
The circles on the bottom faces of the frames indicate the edge of the core. 




Fig. 2. Same fiber and z locations as in Fig.[T]except the amplifier is counter-pumped instead 
of co-pumped. The threshold powers are slightly different for co- and counter-pumped 



pump cladding. We compute the heat deposition rate Q using 

' v„ - V, 



Q{x,y)=NYb{x,y) 



a^,-{a^!, + a'p)nu{x,y) 



Ip, (2) 



where the first term in brackets is the quantum defect. The heat profile can be visualized by 
inverting the depletion profiles in Figs. [T] and |2] and multiplying them by the doping profile 
NYb{x,y). 

The portion of the heat profile that is responsible for mode coupling gain (STRS gain) is the 
antisymmetric part created by the antisymmetric part of the signal irradiance. The population 
saturation strongly influences the shape of this part of the heating. In Figs. [3] and |4] we show 
core centered cuts through the oscillatory heat profiles at the same three z locations as in Figs.[T] 
and 12] These heat profiles create the oscillating asymmetric part of the temperature profile that 
couples modes LPqi and LPn in the STRS process. From an examination of these heat profiles 
it will come as no surprise when we show in the next section that strong saturation causes a 
significant reduction in STRS gain. 

3. Mode coupling gain 

3.1. Without hole burning 

Both Dong and Hansen etal. showed that if the heat profile matches the signal irradiance profile 
over the portion of the core that is (uniformly) doped, the gain of mode LPi i satisfies 

dPix{z) _ 



dz 



n+8oiXPoiiz) Pn{z)^gmtPniz). (3) 



Here we have simplified Dong's expression by assuming there is no depletion of mode LPqi 
due to STRS, and no linear loss for either mode. The total gain gnet for LPn is its laser gain, 
indicated by gi\, plus the STRS gain, indicated by {go\XPoi), where goi is the laser gain for 
LPqi. Thequantumdefectheatdepositedby laser amplification of LPqi is (goiPoi) &o x is a real 
valued coefficient that relates quantum defect heating to STRS gain. For each fiber design x has 
a constant value determined by the frequency offset between LPqi and LPn and the quantum 
defect, plus the thermal, geometrical, and optical properties of the fiber 

Under the assumption that the shape of the heat profile matches the irradiance profile, x is 
independent of z and the modal powers. It depends on the spatial overlap of the two modes with 
one another and with the dopant profile. Hansen et al. and Dong both showed that confined 
doping reduces the value of ;|f , as does reducing the V parameter below 5 or so. 

For full doping and for V > 5, Hansen et al. and Dong both showed that, according to their 
models, if the starting noise level in LPn is set to approximately 10^'^ W corresponding to 
quantum noise, the threshold power lies near 400 W. Details of the fiber design other than V 
and the doping profile do not affect this threshold. Using a starting power of 10^^ W rather than 
10^'^ W reduces the threshold power to approximately 200 W. Confining the doping diameter 
to 50% of the core diameter was found to approximately double the threshold power 

We used our model to verify the value of x for one case analyzed by Hansen et al.{¥\g. 4 of 
ref. El). In order to avoid population saturation we used a pump power of 20 kW, c/core = 40 
/im, fifciad = 250 /im, and an LPqi seed power of 50 W. Otherwise the fiber parameters were 
those from Table 1. Our frequency of maximum gain was equal to that of Hansen et al., and 
our computed a value for x was within 5% of that of Hansen et al. Of course, our fiber was 
extremely inefficient, using 20 kW to produce approximately 400 W of signal, and depleting 
the pump by only 2.5%. However, this exercise does demonstrate good agreement between our 
model and those of Hansen et al. and Dong in the limit of low saturation. 
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Fig. 3. Normalized symmetric part of the heat profile (upper plot) and antisymmetric oscil- 
latory portion of the heat profile (lower plot) for same co-pumped fiber at the same three 
z locations indicated in Fig.[T] Near the input end the symmetric heat profile closely re- 
sembles the LP()i in^adiance profile, while the oscillatory part resembles the product of 
the fields of LPqi and LPn- Farther along the fiber the symmetric part of the heat profile 
becomes nearly flat topped while the antisymmetric part is strongly suppressed near the 
center (x = 0). 
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Fig. 4. Normalized symmetric part of tlie lieat profile (upper plot) and antisymmetric os- 
cillatory portion of the heat profile (lower plot) for the same counter-pumped fiber at the 
same three z locations indicated in Fig.|2] The oscillatory anti symmetric heat profile never 
closely matches the product of the fields of LPqi and LPi i . 



3.2. With hole burning 

The model comparison just described illustrates the lack of realism of the assumption of equal 
shapes for the heat and irradiance profiles. For the heat profile to match the light profile the 
pump irradiance must be much stronger than the signal irradiance. However, amplifiers that 
efficiently convert pump power necessarily experience strong population saturation. To empha- 
size this contrast between our model and one that assumes matching heat and irradiance profiles 
and thus constant X values, we define a new shape factor %' as 



X 



,?comp Ss 
gsPs 



^strs 



(4) 



where gcomp is the total gain of LPi i computed with our model, replacing the gmt of Eq. [3] gs 
is the computed signal laser gain, and gstrs is the STRS gain. We do not distinguish between the 
nearly equal laser gains for the two modes in this high V fiber with fully doped core. Variables 
Ps, gs, and gstre are all z dependent computed values. The denominator is proportional to the 
total deposited heat. Figure|5]shows plots of x' for co-pumped and counter-pumped ampHfiers. 



The upper plot is for a 100 /xm diameter cladding; the lower plot is for a 400 /im diameter 
cladding. Saturation effects are stronger with the larger cladding because of the reduced pump 
irradiances. In both fibers the value of j' near z = for the co-pumped case is nearly equal to 
the % of Dong and of Hansen etal., as expected because of the low degree population saturation 
there. As z increases the population saturation strengthens so the value of j' falls, implying the 
mode coupling gain is reduced relative to the laser gain. This gain reduction raises the STRS 
threshold relative to that of Hansen et al. and Dong. In the small cladding fiber the reduction 
does not occur until half way along the fiber while in the large cladding fiber it occurs much 
sooner 

As expected from the nearly constant ratio of pump to signal for the counter-pumped fibers, 
the value of x' is nearly constant along the fiber However, its value is still reduced to a level 
seen in the co-pumped fiber near the crossing point near z = 0.5 m for the small clad fiber and 
near z = 1 .4 m for the large clad fiber It is clear from a comparison of the values of x' that 
the STRS threshold should be higher for larger cladding sizes, and the threshold for an efficient 
amplifier should exceed the Hansen et al. and Dong threshold. 

Another way to view the same information is to plot §comp, gs, and the total heat rather than 
X'- This is done in Fig. |6]for the same pair of fibers, co-pumped, with rfciad — 100 /zm (upper 
plot) and 400 ^m (lower plot). 

The expected trend of higher thresholds for larger cladding sizes is illustrated in Fig.|2] The 
computed thresholds, defined as 1 % of the signal power in LPj i , are seen to rise with increas- 
ing cladding diameters and the resulting increasing degree of population saturation. Co- and 
counter-pumped fibers are found to have similar thresholds. Decreasing the doping diameter 
also increases the threshold because saturation is then strong across the full doping profile. For 
comparison we show as a black horizontal line in the figure the threshold predicted without 
saturation. More details for the model runs included in Fig. |2]are listed in Tables 2-5. 

4. Stimulated Brillouin scattering (SBS) suppression 

Designing a narrow bandwidth fiber amplifier suitable for beam combining applications re- 
quires balancing SBS and STRS. Here we present a simplified version of this balancing act. We 
define an effective SBS gain length at the STRS threshold by 

^thres-Leff= / Ps{z)dz. (5) 

Jo 
The SBS threshold power |[TOl is exceeded if 

gB fthres^eff 
7 ^eff 



> 17, (6) 



where the 17 comes from the usual SBS threshold gain value of 21 minus the laser gain for the 
Stokes wave ifTOl . The 7 parameter is an SBS gain reduction factor that reduces the effective 
SBS gain from its nominal value of ^b = 5 x 10^" mAV for silica. The linewidth for silica is 
approximately 50 MHz and the Stokes shift is approximately 16 GHz for X = 1032 nm. The 
value of 7 can be increased by broadening the signal linewidth by phase modulation jToll above 
50 MHz, for example, or by introducing a temperature gradient 1 11 1 or a strain gradient I112II131 
along the fiber length to vary the Stokes shift. Another approach is to vary the acoustic velocity, 
and thus the Brillouin shift, across the core region ||T4][T51 . Temperature and strain shift the 
SBS frequency by known amounts, so from the signal linewidth and the temperature and strain 
gradient one can estimate 7 using an SBS model. To avoid SBS the value of 7 must satisfy 
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Fig. 5. Plots of "X^ versus z for co-pumped (dashed green curve) and counter-pumped (solid 
blue curve) fibers operating near the mode instability threshold. The fiber parameters: 50 
/im diameter core and doping, 100 \lra diameter pump cladding (upper plot) and 400 ^m 
diameter pump cladding (lower plot), 11 00 Hz red detuning of LP; i , A^ = 1032 nm, Xp = 
976 nm, A'A=0.054 («core = 1-451, ^clad = 1-45). In the upper plot the pump powers are 525 
W co-pumped, and 493 W counter-pumped. In the lower plot the pump powers are 1200 
W co-pumped, and 1350 W counter-pumped. 
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Fig. 6. Plots of laser gain gj, total gain gcomp and heat versus z for the same co-pumped 
fibers and the same operating conditions as in Fig. |5] The upper plot is for a 100 jlva 
diameter cladding; the lower plot is for a 400 ^m diameter cladding. In both plots it is 
clear that the gain and heat profiles are not closely matched. 



For our 50 /im diameter core, Aeff = 1 175 /im^, so the SBS threshold condition is 



7 = 



Mhres^eff 

400W-m' 



(8) 



If Pthres^efT > 400 W-m, this expression gives the minimum value of y necessary to suppress 
SBS. 

In Fig. [8]we plot the STRS threshold powers versus the quantity Pthies^efT for our example 
fiber with c/core = 50 /im. Selecting an arbitrary value for Ahresi'eff on the horizontal axis (or 
equivalently the value of 7), the values of the STRS threshold can be read from the computed 
curves. From the four curves it appears that for any value of 7 the counter-pumped fully doped 
amplifier provides the highest STRS threshold power, the co-pumped fully doped amplifier 
offers the lowest STRS threshold power, while the confined doping, co-pumped amplifiers give 
intermediate threshold powers. 



5. Scaling for other core/cladding sizes 

To first order changing the core diameter while keeping the ratio t/core/'^ciad fixed, and keeping 
NA = 0.054, the threshold powers do not change. The frequency offset changes proportional to 
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Fig. 7. Output signal powers at the instability threshold for different pump cladding diame- 
ters (100, 200, 300, 400, 500 ^um). The parameters: dcore = 50 /xm, X,, = 976 nm, Aj = 1032 
nm, signal seed powers are 10 W in LPgi and 10"'^ W in LPn. The lengths of the fibers 
are adjusted the minimum value necessary to achieve high efficiency, defined as pump ab- 
sorption > 0.95. The solid curve at 345 W indicates the threshold computed in the limit of 
no population saturation. Details are given in the tables of Section[8] 
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Fig. 8. Plot of signal threshold output power versus -Pttire-s^eff for the 50 flm core fiber. The 
curves are spline fits to the computed data points. For Pthres^eff = 2400 W-m (7 = 6) the 
threshold powers are 1050 W for co-pumped and fully doped; 1205 W for co-pumped and 
^^dope = 40 flm; 1347 W for co-pumped and d^ope = 30 flm; 1572 W for counter-pumped 
and fully doped. 



1/Aeff but the degree of saturation, and thus the value of j' is unaltered. 

However, for Jcore < 30 /im, the value of Y is reduced to less than 5, so the mode overlap 
with the core becomes noticeably weaker, which tends to raise the STRS threshold ||6ll8l. At the 
other extreme, for dcore > 80 /im, thermal lensing becomes significant, and this constriction of 
the modes may alter the STRS threshold. We will examine these cases more closely in future 
studies. 

As rfcore changes the area Aeff changes and this changes the value of 7 required to suppress 
SBS for a given value of Pthres^eff according to Eq. |8] Similar powers are possible, while avoid- 
ing both SBS and STRS, but the value of 7 necessary to avoid SBS scales approximately as 
V'^core- ^f there is a constraint on 7 imposed by beam combining requirements, for example, 
this implies a lower limit on the core size. 

6. Scaling for other LPi 1 starting powers 

The thresholds computed in this study are based on an input signal power of 10^'^ W in LPn, 
with a frequency shift to the STRS gain maximum, or threshold minimum. This input power 
level corresponds approximately to the quantum noise limit, and serves as a standard point 
of comparison among models. Actual starting powers will be considerably higher than this if 
there is amplitude modulation on the pump or signal seed with a modulation frequency near the 
optimum STRS frequency shift JSjIll. Thermal noise is probably 2-3 factors of 10 higher than 
this as well. More realistic thresholds would be based on input powers of perhaps 10^"* W, but 
the exact level will require measurements of the modulation properties of the pump and seed. 
The reference thresholds presented here (Pref) can be used to estimate the threshold for a higher 
input power (Pstart) using 

log(fitart/10) 

'log(lO-iVlO) 

where the divisors of 10 are present because our definition of threshold as 1% of the power in 
LPn implies a threshold output power of order 10 W in LPn . 

The thresholds computed in this study are also based on negligible signal loss due to absorb- 
ing impurities in the glass or due to photodarkening. The presence of these processes adds heat 
with a profile that contributes to STRS gain, and so reduces the thresholds 0. 

7. Conclusion 

We showed that transverse hole burning, or saturation of the population inversion, in a Yb^+ 
doped fiber can strongly influence the STRS or mode coupling gain, and this effect can be 
exploited to substantially raise mode instability thresholds. Saturation can also be used in con- 
junction with confined doping, both contributing to raising the threshold. 

The benefits of high saturation persist when SBS suppression is considered, even though high 
saturation implies longer fibers with higher Pthresi-eff values. The benefits of confined doping 
also survive the requirement of SBS suppression. 

Fibers with high saturation have the added advantage of minimizing photodarkening, assum- 
ing photodarkening increases with higher upper state population density, as experiments seem 
to indicate lfT6l . 



thres — Pref,__/-,r,-16 /1/AN '"' 



8. Tables 

In the following tables Av is the frequency of LPqi minus the frequency of LPi i . 

Table 2. Thresholds: co-pumped, dcore = 50 ^m, ddope = 50 M™ 
«fciad[Mm] Av[Hz] L[m] Pthres [W] JPdz[W-m] L,s 



100 


1100 


0.8 


488 


161 


0.330 


1 


200 


1100 


1.6 


685 


603 


0.880 


1.51 


300 


1100 


2.6 


885 


1375 


1.55 


3.44 


400 


1100 


4.0 


1101 


2789 


2.53 


6.97 


500 


1100 


6.0 


1335 


5338 


4.00 


13.3 


Table 3 


. Thresholds 


: counter 


-pumped, dcore 


= 50 /xm, rfdope 


= 50 iim 




'^clad [Mm] 


Av [Hz] 


L[m] 


Ahres [W] 


JPdz [W-m] 


Leff 


7 


100 


1100 


0.8 


453 


98 


0.216 


1 


200 


1100 


1.6 


676 


281 


0.416 


1 


300 


1100 


2.6 


921 


634 


0.688 


1.58 


400 


1100 


4.0 


1220 


1288 


1.06 


3.22 


500 


1100 


6.0 


1580 


2429 


1.54 


6.07 



Table 4. Thresholds: co-pumped, d^ore = 50 flm, rfjope = 40 flm 
Jdad[Aim] Av[Hz] L[m] P.hres [W] JPdz[W-m] L.g 



100 


1400 


0.9 


549 


204 


0.372 


1 


200 


1400 


1.9 


844 


893 


1.06 


2.23 


300 


1400 


3.2 


1185 


2290 


1.93 


5.72 


400 


1400 


5.0 


1567 


4974 


3.17 


12.4 



Table 5. Thresholds: co-pumped, dcore = 50 /im, ddope = 30 flm 
Jciad[Mm] Av[Hz] L[m] Ahres [W] JPdz[W-m] Us 



100 


1900 


1.4 


786 


493 


0.627 


1.23 


200 


1900 


3.0 


1311 


2257 


1.72 


5.64 


300 


1900 


6.0 


1975 


7757 


3.93 


19.4 



